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Abstract 

The global attractor conjecture says that toric dynamical systems (i.e., 
a class of polynomial dynamical systems on the positive orthant) have 
a globally attracting point within each positive linear invariant sub¬ 
space - or, equivalently, complex balanced mass-action systems have 
a globally attracting point within each positive stoichiometric com¬ 
patibility class. A proof of this conjecture implies that a large class of 
nonlinear dynamical systems on the positive orthant have very simple 
and stable dynamics. The conjecture originates from the 1972 break¬ 
through work by Fritz Horn and Roy Jackson, and was formulated in 
its current form by Horn in 1974. 

We introduce toric differential inclusions, and we show that each 
positive solution of a toric differential inclusion is contained in an in¬ 
variant region that prevents it from approaching the origin. We use 
this result to prove the global attractor conjecture. In particular, it fol¬ 
lows that all detailed balanced mass action systems and all deficiency 
zero weakly reversible networks have the global attractor property. 
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1 Introduction 


Any autonomous polynomial dynamical systemQon the strictly positive or- 
thant M” can be represented as 

- y) ( 1 ) 

y^y'&G 


for some geometrically embedded graplij^ G and some positive constants 
ky^yi, one for each edg^ y ^ y' of G. (Here denotes the monomial 
Similarly, any non-autonomous polynomial dynamical system on 
can be represented as 

h^y'it)x^iy' -y) (2) 

y^y'eG 

for some nonnegative scalar functions ky^y'{t). 

Positive trajectories of ([^ and ([^ are also trajectories of mass-action 
systems na. There is great interest in understanding the persistence and 
global stability properties of such dynamical systems [ni[i5i[9i[ini[i2i[iii[i9i 
12niEni[2Ill221ESl[niEl[21EHll3ll3[ni[Il0[in]. a natural question is the 
following: for which systems ([^ or (|^ is it true that trajectories that start 
in stay awajj^from the boundary of M”? 

Given a polynomial dynamical system, the choice of graph G above is 
not unique. We will see that if the graph G can be chosen such that each 
edge of G is contained in a cycl^ then any bounded trajectory of ([^ must 
be persistent in M”. More generally, this is also true for solutions of (|^, 

^or more generally any power law system [T], since we do not assume that the coor¬ 
dinates of the vertex points y are either integer or non-negative. In particular, if vertex 
points have non-negative integer coordinates we obtain mass-action systems. 

^ A geometrically embedded graph G is a finite directed graph whose set of vertices is 
a finite set Y C K", and each edge of G is represented by an oriented line segment that 
connects two vertices y,y' € Y. See Fig. for examples. 

^Inspired by notation from reaction networks, we denote an oriented edge from y to y' 
by y —> yb Moreover, if y —>■ y' is an edge in G we simply write y ^ y' G G. 

^ i.e., for a trajectory x{t), we have liminf(I) > 0 for all i = Then we say 

i—>-oo 

that the trajectory x(t) is persistent. 

^ or, equivalently, each connected component of G is strongly connected. If this is the 
case, we will say that the graph G is weakly reversible. 
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Figure 1: Four examples of geometrically embedded graphs in Note that, 
while we do assume that the points yi,...,yk G are distinct, we do not 
assume that the line segments (i.e., arrows) representing the vectors yj — yi 
are disjoint. The two graphs on the right are weakly reversible, and the 
bottom-right graph is reversible. Although the bottom-left graph is not weakly 
reversible, it generates dynamical systems which can also be represented 
by some weakly reversible graph. This is not true for the top-left graph. 
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provided that the nonnegative scalar functions ky^y/{t) are bounded away 
from zero and infinit'^ 

In this paper we prove the results mentioned above, and use them to 
prove the global attractor conjecture, which says that toric dynamical sys¬ 
tems have a global attractor within any linear invariant subspac^ A toric 
dynamical system is a polynomial dynamical system|^ for which there exists 
a representation of the form Q that admits a vertex-balanced equilibrium in 
M”, i.e., there exists xq G Wf such that for each vertex y eY we have 


E u _^y — h- 

rVy^yXQ — / ^ IVy^yXQ. 

y^y&G y^y&G 


(3) 


In other words, if we think of the positive number ky^yix^ as a flow from the 
vertex y to the vertex y', then condition ([^ says that at each vertex of the 
graph G, the sum of all the incoming flows equals the sum of all outgoing 
flows. The identity ([^ implies that the graph G is weakly reversible [IT] . 

The name “toric dynamical system” has been introduced recently |H] to 
emphasize the remarkable algebraic properties of these systems, but this 
class of polynomial dynamical systems has been hrst studied in depth in the 
1972 breakthrough paper of Horn and Jackson la, where they have been 
called complex balanced mass-action systems. Also in [IT] Horn and Jackson 
have shown that these systems enjoy remarkable stability properties, and in 
particular they have a unique positive equilibrium within each linear invari¬ 
ant subspace, and this equilibrium is locally asymptotically stable. Actually, 
in na Horn and Jackson stated that the unique positive equilibrium within 
each linear invariant subspace is a global attracto^ but soon afterwards Horn 
explained that they have not actually proved this claim, and in 1974 he pro¬ 
posed this global convergence property as a conjecture [T3] . 

The conjecture that toric dynamical systems have a global attractor within 
each linear invariant subspace (or, in the language of Horn, that complex 
balanced mass-action systems have a global attractor within each reaction 
simplex) was later called the Global Attractor Conjecture [8]. 

® i.e., there exists £ > 0 such that ky^yft) S [e, j] for all t. 

^ A linear invariant subspace of Q is an invariant set Sq = {xq -+- Sq) H M" , where Sq 
is a linear subset of K”. For toric dynamical systems we will see that S'o = span{y' — y : 
y ^ y' £ G} gives rise to a linear invariant subspace for any xq £ Rif- 

®or more generally a power law system 

® i.e., if the positive equilibrium xq belongs to the linear invariant subspace So, then 
all trajectories that start in Sq converge to xq- 
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The global attractor conjecture has resisted efforts for a proof for over 
four decades, but proofs of many special cases have been obtained during 
this time, for example [HI EHl El El El EZl HI IT6] . 

In particular, Craciun, Nazarov and Pantea [1] have recently proved the 
three-dimensional case, and Pantea has generalized this result for the case 
where the dimension of the linear invariant subspaces is at most three [TB] . 
Using a different approach, Anderson has proved the conjecture under the 
additional hypothesis that the graph G has a single connected component , 
and this result has been generalized by Gopalkrishnan, Miller, and Shin for 
the case where the graph G is strongly endotactic [TB] . 

The results described above do not provide a proof of the global attractor 
conjecture in full generality, and also do not provide a proof for any of the 
following three important special cases: [i) the case where G is reversiblG^ ^ 
(ii) the case where G is reversible and the system Q is detailed balanceG^, 
and {Hi) the case where G is reversible and ky^yi = 1 for all y ^ y' E G. 

The goal of this paper is to introduce toric differential inclusions, and 
to use them to construct a proof of the global attractor conjecture in full 
generality. 


Toric differential inclusions are defined as follows. Consider a finite set 
of polyhedral cones^ that cover M”, such that is a polyhedral faij^ 
[271 EH]. For each cone G E T denote by G° the polar cona^of G. Fix some 
5 > 0. For each x E MT define Fjrs{x) to be the convex cone generated by 
the union of polar cones G° for aX\ G E F such that dist{x, G) < 6 . Then a 


i.e., ii y ^ y' E G then also y' ^ y E G. 

i.e., we replace the vertex-balance equilibrium condition ([^ with the more restrictive 
edge-balance equilibrium condition: ky^yXQ = ky^yxfj for all y ^ y E G. This condition 
says that, at the equilibrium point xq, the “forward flux” balances the “reverse flux” for 
any edge in G. 

A polyhedral cone in M” is the intersection of a finite set of half-spaces of K". 

A polyhedral fan in M” is a finite set T of polyhedral cones such that (i) any face 
of a cone in T is also in T, and (ii) the intersection of two cones in is a face of both 
cones. We say that a polyhedral fan T covers K” if Simple examples of 

polyhedral fans Jhr in K" are given by the polyhedral cones delimited by a finite set H of 
hyperplanes through the origin. This particular case of “hyperplane-generated polyhedral 
fan” is especially relevant for motivating our definition of toric differential inclusions. 

The polar cone of a cone G E K" is the cone G° = {y E K"| x-y < 0 for all x E G} |29) . 
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toric differential inclusion is a differential inclusion on Wf given by 

dx ^ , 

for some polyhedral fan iF as above, and some 5 > 0. 

It turns to that, although the global attractor conjecture is formulated in 
terms of systems if the form Q, we need to understand some properties of 
systems of the form (|^ in order to construct a proof of the conjecture. There¬ 
fore, we hrst prove that we can reduce some systems (|^ to toric differential 
inclusions: 

Theorem A. If a polynomial dynamical system has a representation ^ 
such that G is weakly reversible and the non-negative functions ky^yft) are 
bounded away from zero and infinity, then it can be embedded^ in a toric 
differential inclusion. 

In particular, it follows that if G is weakly reversible, then Q can be 
embedded in a toric differential inclusion for any choice of parameters ky^y' > 
0. We then prove that, for any toric differential inclusion T and for any 
positive point xq € M”, there exists a hypersurface H C that separates 
into two regions, such that 0 and Xq are in different regions, and the region 
that contains the point Xq is an invariant region for the solutions of T: 
Theorem B. Toric differential inclusions have an exhaustive set of zero- 
separating hypersurfaces. 

Then we use this result to conclude that positive trajectories of toric 
dynamical systems must stay away from the boundary of : 

Theorem C. Toric dynamical systems are persistent. 

Finally, we are able to conclude that the global attractor conjecture is 
true: 

Theorem D. Toric dynamical systems have a globally attracting point within 
each positive linear invariant subspace. 

2 Definitions and notation 

As we mentioned in the previous section, some of the main objects of in¬ 
terest in this paper are non-autonomous polynomial dynamical systems of 

We say that the dynamical system ^ = f{x) is embeddedinto the differential inclusion 
G Fix) in the domain O, if f{x) G F{x) for all x G ti. 
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the form ([^ such that the graph G is weakly reversible, and the nonnega¬ 
tive scalar functions ky^y>{t) are bounded away from zero and inhnity. We 
will refer to this class of dynamical systems as k-variable toric dynamical 
system^^ 

More precisely, a k-variahle toric dynamical system is a dynamical system 
on of the form 


-y) ( 4 ) 

y^y'&G 


such that G is weakly reversible, and there exists a hxed positive number e 
with e < ky^y'it) < 4 for all t and for all y ^ y' E G. 

Note that, under these conditions, solutions of the fc-variable toric dy¬ 
namical system above are also solutions of the differential inclusion on M” 
given by 


— G F{x), where F{x) 


ky^y,xy{y'-y) 

y^y'^G 


E ^ ky^yf 



( 8 ) 

We will show that solutions of the differential inclusion (|^ are also solu¬ 
tions of a type of differential inclusions with rich geometric properties, called 
toric differential inclusions. Later in this section we will dehne the general 
notion of toric differential inclusions in M”. We hrst dehne the class of polar 
differential inclusions in M”. 


2.1 Polar differential inclusions 


If the graph G is reversible, then the equations Q can be written as 


dx 

dt 


Y iy' - y)^ 

y^y'eG 


( 6 ) 


by grouping together terms given by an edge y ^ y' and its reverse y' ^ y. 
Then, for each reversible edge y ^ y' E G, there is now a single term 

The name k-variable toric dynamical system is due to the fact that any such dynamical 
system can be obtained from a toric dynamical system by using the same graph G, and 
replacing the constants ky^yi with time-dependent non-negative functions ky^y/ft) that 
are bounded away from zero and infinity. 
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in the sum (§, which can be thought of as a “tug-of-war” between the 
forward and reverse terms. Indeed, both the forward and the reverse terms 
are trying to “push” the state x{t) of the system along the same lin^^ but 
in opposite directions. Then, the domain can be partitioned into three 
regions: the region where the inequality ex^ > ^x^ holds (which implies 
ky^y/{t)xy > ky'^y{t)x^'), the regiou where the inequality ^x'^ < ex'^' holds 
(which implies ky^yi{t)x^ < kyi^y{t)x'^'), and an uncertainty region where 
neither one of these two inequalities are satished, and either one of the two 
terms ky^yi{t)x^ and ky'^y{t)x^ may win the tug-of-war, due to the fact that 
ky^yi{t) and ky/^yit) may take any values between e and -. Some concrete 
examples can be found in Section 3 of |1] . In particular, note that if one term 
does win the tug-of-war, then the direction of the sum of the two terms is 
towards the uncertainty region, and also towards the hypersurface x'^ = x"^ 
in . 

The simplest way to understand how these regions adjoin each other is 
to look at them in instead of M”, after applying a logarithmic transfor¬ 
mation. If we denote X = log x, then the hypersurface x"^ = x"^ becomes the 
hyperplane {y' — y)-X = 0. Moreover, the uncertainty region described above 
becomes just the set of points at distance < 6 = from this hyperplane. 

Note also that the direction given by the vector y' — y is exactly the direction 
orthogonal to the hyperplane {y' — y) ■ X = 0. 

Therefore, if the graph G consists of a single reversible edge y y', then 
the dynamics of the system ([^ at a point x can be described as follows: we 
know that x will move along a line of support vector {y' — y), and we are able 
to specify one direction or another by mapping x to X = log x, and checking 
whether the distance between X and the hyperplane {y' — y) ■ X = 0 is > S. 
Moreover, if we can specify this direction, then it is always the direction 
towards (and not away from) this hyperplane, and orthogonal to it. 

This characterization has the advantage that it can be carried over to the 
more general case, where G contains several reversible edges. In that case we 
have several tug-of-wars going on at the same time, but for each one of them 
we can specify the winning direction (if any) at x by calculating the distance 
between X = logx and some hyperplane in M"’. Depending on whether X 
falls outside an uncertainty region or not, each reversible edge y ^ y' oi G 
contributes one or two resultant vectors (if one, then it is either y' — y or 
y — y', and if two, then they are ±{y' — y)). 


i.e., a line whose direction is given by the vector y' — y. 




If follows that any system ([^ can be embedded into a differential inclusion 
on given by a set Ti of hyperplanes in M"' and a number 5 > 0, as follows. 
For each x G we dehne F-^i^s (Jog x) to be the convex cone generated by 
vectors orthogonal to the hyperplanes of 'H, in the direction that goes from 
the point X = logx towards each hyperplane, and also the opposite direction 
if X is at distance < 6 from some hyperplane. If X does not belong to any 
uncertainty region, then F-^^siJogx) is dehned to be exactly the polar cone 
C° of a cone C bounded by a subset of hyperplanef^ from "H. Moreover, if X 
does belong to some uncertainty regions, then we can still describe F-^, 5 (log x) 
in terms of polar cones, by including not just the polar of the largest cone 
of F{'H) that contains X, but also the polar of each cone of F{'H) that is at 
distance < 5 from X. 

Of course, not all polyhedral fans are determined by a set of hyperplanes 
as above. Nevertheless, we can generalize the construction described above 
to dehne polar differential inclusions given by a polyhedral fan F in M"’, 
as described in the previous section: we dehne Fjrs(X) to be the convex 
cone generated by the union of polar cones C° for all O' G X such that 
dist{X, C) < 6. 

2.2 Toric differential inclusions 

Polar diherential inclusions can produce sets of cones that contain the right- 
hand-side of the vector helds Q for reversible G, but they are dehned on M"’, 
while the system Q is dehned on . To obtain a proper generalization we 
need to introduce toric differential inclusions, which are obtained from polar 
diherential inclusions by a logarithmic change of variable: a toric diherential 
inclusion given by a polyhedral fan F in is a diherential inclusion on M” 
given by 

dx ^ 

e F^^si^ogx) (7) 

where Fj^ s is a polar diherential inclusion as dehned above. 

Due to the way this dehnition is related to the system ([^, it follows that 
fc-variable toric dynamical systems (|^ can be embedded into polar diherential 
inclusions, provided that the graph G is reversible. 

the cone C is the largest cone that contains the point X within the polyhedral fan 
F{7i) determined by the set of hyperplanes H. The convex cone generated by the “outer 
normal” vectors orthogonal to the hyperplane faces of C is its polar cone C° [55] • 
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As we discussed above, for a reversible graph G, the polyhedral fan 
could be chosen simply as the fan generated by the set of hyperplanes that 
are orthogonal to the vectors y' — y ioi y ^ y' ^ G. 

We need to remove this strong reversibility assumption, in order to be 
able to address the global attractor conjecture in full generality. In the next 
section we explain how to hnd such a polyhedral fan for any weakly reversible 
graph G. 

3 Connection between /c-variable toric dynam¬ 
ical systems and toric differential inclusions 

As we discussed in the previous section, the simplest examples of polar differ¬ 
ential inclusions are generated by polyhedral fans that are determined by 
a hnite set T-L of hyperplanef^ that pass through the origin. Then, the sim¬ 
plest examples of toric differential inclusions are also generated by polyhedral 
fans J-y, as above. We will refer to this class of toric differential inclusions as 
hyperplane-generated toric differential inclusions. 

We have seen in the previous section that any reversible /c-variable toric 
dynamical system in can be embedded into a (hyperplane-generated) 
toric differential inclusion. Here we show that a similar fact is true for all 
fc-variable toric dynamical system^^ 

Theorem 3.1. Consider a k-variable toric dynamical system ^ . Then this 
system can be embedded into a toric differential inclusion. 

Proof. Assume hrst that the weakly reversible graph G is made up of a single 
oriented cycle. If our single-cycle graph is given by j/i —)■ 2/2 —t ... —t |/r —t |/i, 
then the /c-variable toric dynamical system it generates is of the form 

dr 

= ^ki{t)xy^{yi+i-yi), ( 8 ) 

2=1 

where yr+i = yi and Eq < kiif) < — for some £0 > 0. 

i.e., Py consists of all the cones delimited by hyperplanes in 71. 

In order to prove the global attractor conjecture we need to construct zero-separating 
surfaces for all fc-variable toric dynamical systems (not just for the reversible ones). Theo¬ 
rem implies that it is sufficient to construct zero-separating surfaces for toric differential 
inclusions. 
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Consider the set C of lines throngh the origin in the direction of vectors 
Hi — Uj for all i ^ j, and denote by Ti the set of all hyperplanes that are 
orthogonal complements of lines in Denote by the polyhedral fan 
generated by the set of hyperplanes Ti, and, for 5 > 0 , denote by Tn,s the 
corresponding hyperplane-generated toric differential inclnsion. We will show 
that there exists 5o > 0 snch that the single-cycle /c-variable toric dynamical 
system (|^ is embedded in the toric differential inclnsion Th, 5 o- 

Choose ^0 > 0 large enongh snch that the nncertainty regions given by 
the reversible edges yt ^ yj and Eq are contained within the nncertainty 
regions of the toric differential inclnsion Tu^So- 

Fix a point x G M" snch that X = logx belongs to some cone C in 
Tu-, and does not belong to any nncertainty region of 7 ^, 50 • particnlar, it 
follows that the cone C has dimension n, otherwise X wonld be contained in 
some hyperplane in which in tnrn wonld be contained in an nncertainty 
region of T'n, 5 o- We want to show that the right-hand-side of (|^ is contained 
in the dnal cone C°. 

Consider a vector w in the interior of C, and project yi,y 2 , ■■■,yr on the 
line Iw that passes throngh the origin in the direction given by w. Then no two 
projections are the same, becanse w does not belong to any of the hyperplanes 
in l-i. We now give a second set of names to the vectors yi,y 2 , ■■■,yr, say 
Vi,V 2 , ...,Vr, to record the distance to the origin of the projections along the 
line (i.e., if the projection of yi^ is the closest to origin, then vi = and if 
the projection of yi^ is the second closest to origin, then V 2 = yi 2 , and so on). 
Note that, since the interior of C is disjoint from all the hyperplanes in T-L, it 
follows that the new names ni, ^ 2 , Wr do not depend on the particnlar choice 
of vector w in the interior of C. In other words, the dot prodncts {vi+i—vi)-w 
are all negative nnmbers, for all w in the interior of C. Therefore, the vectors 
V 2 — ni. Vs — V 2 , ...,Vr — Vr -1 belong to C°. 

So, in order to show that the right-hand side of ([^ is inclnded in C°, it 
is enongh to show that it can be written as a positive linear combination of 
the vectors V 2 — fi, ^3 — V 2 , •••, Vr — Wr-l- 

Note that (ii, * 2 , •••, v) is a permntation of (1,2, If we denote the 

inverse permntation by ••■,>), h follows that yi = y 2 = and 

so on. 

a hyperplane belongs to % iff it is orthogonal to a line in C. 
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Then we have y2 — Vi = Vj^ — Vj^. If j2 > ji we write 

j2-i 

2/2 - 1/1 = ^{vi+i-vi), 

l=h 

and if j 2 < ji we write 


ii-i 

2/2 - 2/1 = - -vi), 

l=h 

We do the same for 2/3 — 2 / 2 , 2/4 — 2 / 3 , and so on. This way, we write each 
difference 2/i+i — Vi from the right-hand-side of ([^ in terms of the vectors 
±( 2 /i+i — 2/z), with I = 1,2, ...,r — 1. Therefore we can re-group terms to 
obtain ^ 

^ (9) 

i=l 

where is a sum of several terms of the form ki with various signs. 

Note now that the positive terms inside <I>; correspond to edges of the form 
Vm — t Vn with m < I < n, and negative terms inside correspond to edges 
of the form Vm —t Vn with n < I < m. This means that the positive terms 
inside contain kix'^^ with i < I, and the negative terms inside <I>z contain 
kix'"' with i > 1. Since logo: G C (and is not in an uncertainty region), and 
due to our choice of ^o, if follows that kix^^ < k 2 x'^^ < ... < kj.x'"^. Therefore, 
the sum of the positive terms inside dominates the sum of the negative 
terms inside <h;, for each Q In conclusion, the right-hand-side of ([^ (and 
therefore ([^) is a positive linear combination of the vectors vi+i — vi, for 
1 < Z < r — 1, so it belongs to C°. 


If X does belong to an uncertainty region of the toric differential inclusion 
T'h,So^ then denote by F{x) the cone of Th, 5 o at It follows that the maximal 
linear subspace contained in F{x) has dimension > 1 in M”. We project the 
problem on the orthogonal complement of that subspace, and then we reason 
the same way as abovt^^ This gives us the desired conclusion for the case 


when G is made up of a single oriented cycle. 


note that the number of positive terms inside is the same as the number of negative 
terms inside because the graph G is a cycle. 

When projecting on the smaller dimensional subspace, the graph G is replaced by its 
projection G. All vertices of G that project to the same vertex of G are interchangeable 
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If the weakly reversible graph G is not a single oriented cycle, then we 

9 

write it as a union of cyclic graphs, G = and we can argue as above 

i=l 

for each such Gi. We obtain that the /c-variable toric dynamical systems 
given by the cycle Gi are embedded in toric differential inclusions generated 
by some set of hyperplanes Hi. Note now that the right-hand-side of a k- 
variable toric dynamical system given by G can be decomposed into a sum of 
terms, such that each term is of the formp^ given by the right-hand-side of a 
fc-variable toric dynamical system determined by G*. Then we conclude that 
any fc-variable mass-action system given by G can be embedded into a toric 

9 

differential inclusion generated by the set of hyperplanes H = [^Hi. □ 

i=l 

4 Construction of zero-separating surfaces for 
toric differential inclusions 

In this section we show that for any toric differential inclusion T on M” 
and any small enough neighborhood Vq of the origin in M"" there exists a 
zero-separating surface, i.e., a hypersurface ^r,Vo such that \ is the 
union of two disjoint and connected open sets Z^y^ and Z^y^, such that 
Z^y^ C Vo, Z\- is an invariant region of T, and the closure of Z\- does not 
contain the origin. 

Then, it will follow that any solution of T with initial condition xq G M” 
must be contained in some invariant region Zf- as described above, and in 
particular it cannot have an cn-limit point at the origin. 

Remark 4.1. Note that, if we want to construct zero-separating surfaces 
for toric differential inclusions in M", it is actually sufficient to focus on 
hyperplane-generated toric differential inclusions, because for any toric dif¬ 
ferential inclusion Ti there exists a hyperplane-generated toric differential 
inclusion T 2 , such that Ti is embedded in T 2 (for example, we can choose T 2 
to be generated by the set of hyperplanes H that consists of all the support 

with each other when checking that the right-hand-side of ([^ is contained in the (degen¬ 
erate) cone F(x), because the projections are done along linear subspaces contained in 
F{x). 

We may have to use smaller values for the terms in the decomposition, because 
the same edge of G may belong to several graphs Gi. 
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hyperplanes of the {n — 1)-dimensional faces of the cones ofTi). Therefore, 
any zero-separating surface for Ti will also be a zero-separating surface for 

Ti. 

Remark 4.2. Similarly, it is also sufficient to focus on toric differential 
inclusions generated by a set of hyperplanes Ti such that the span{C) = M”, 
where C is the set of orthogonal lines for the hyperplanes in Ti. Indeed, if 
a toric differential inclusion T has span{C) 7 ^ MT, then we can add to the 
set Ti one or more hyperplanes to obtain Ti such that span{C) = Now 
notice that the cone F{x) of the new toric differential inclusion T contains 
the cone F[x) of the old toric differential inclusion T, for all x G M". The 
inclusion relationship between the cones is due to the fact that all generators 
of a cone of T are contained among the generators of a corresponding cone 
of T (actually, the old cones are exactly the intersections between the new 
cones and the space span{C)). Therefore T is embedded in T. 

4.1 One-dimensional toric differential inclusions 

Consider a toric differential inclusion denoted 7i in IR+. Then, close enough 
to the origin, 7i is simply the differential inclusion ^ > 0. In this case a 
zero-separating surface is just a point Pq G M+ such that the cone of 7i at 
Pq points away from zero. But, any point in IR+ that is close enough to zero 
has this property. 

4.2 Two-dimensional toric differential inclusions 

Consider a toric differential inclusion denoted F 2 in 1^+- Then, without loss 
of generality we can assume that F 2 is determined by a set of half-lines that 
start at the origin and bound cones in that can be used to generate 72 . 
More exactly, for each such cone C its exponential image expiC) is a subsefp^ 
of where F 2 is constant and equal to the polar cone C°, except near the 
boundary of exp{C), where, far enough from the point (1,1) 72 is a half¬ 
space, and near the point (1,1) 72 is the whole space Note that (as we 
can see in Fig. in the construction of a zero-separating surface for 72 only 
the half-lines contained in the 3rd quadrant of will really play a role. 

There are many different ways to construct a zero-separating surface 
(which in this case is a zero-separating curve) for 72 . The simplest way 

we will refer to such as set exp{C) as an exp-cone of Tz- 
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Figure 2: A polyhedral fan in (a) gives rise to a toric differential in¬ 
clusion in (b). Neighborhoods of the curves shown in (b) delimit the 
uncertainty regions of the toric differential inclusion. We do not see these 
neighborhoods in (b), but we should imagine that each curve in (b) has some 
nonzero thickness that represents its uncertainty region, and the cone of the 
toric differential inclusion within that uncertainty region is a half-plane (each 
such half-plane is the polar cone of a half-line in (a)). This imposes that the 
zero-separating curve shown in (c) crosses these curves along line segments 
of specified slope; for each crossing of a curve in (c) the line segment must 
be orthogonal to the corresponding half-line in the 3rd guadrant of (a). In 
(d) we see that we may change the slope of the zero-separating curve outside 
the uncertainty regions. This is needed in order to obtain a faithful zero- 
separating curve. See also Section 3 in for related examples discussed in 
more detail. 
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is to construct a polygonal line that connects a point on the x-axis to a 
point on the y-axis, in such a way that there exists exactly one polygonal 
vertex within each bounded exp-cone of 72, and the direction of each line seg¬ 
ment follows the attracting direction of the (single) uncertainly region that 
it intersects (see Fig. 

The construction described above is essentially the same as the construc¬ 
tion of an invariant region in [1], but without having to continue the polygonal 
line all around the point (1,1) to obtain an invariant polygon. 

On the other hand, note that we could have built the polygonal line 
above by using less stringent requirements on its line segments. For example, 
it is not necessary that each line segment follows the attracting direction 
of some uncertainly region; as long as this property is satished within that 
uncertainty region, then in all other regions we can have several line segments 
in various directions, as long as each one of them is a support line of the 
(locally constant) cone of T 2 in that region. Moreover, the zero-separating 
curve does not have to be a polygonal line; it can also contain arcs with 
nonzero curvature, as long as they are connected by straight line segments 
when crossing uncertainty regions (see Fig. j^d)). 

Note also that we can assume that the zero-separating curve constructed 
as above has the following additional property: within each bounded exp-cone 
of 72, the slope of the tangent line to the curve at points outside uncertainty 
regions belongs to the interior of the interval bounded by the slopes of the 
attracting directions of the uncertainty regions. We will say that such a curve 
as a faithful zero-separating curve. This special type of zero-separating curve 
will be used in the next section. 

Note that since we can prove, as above, that two-dimensional toric dif¬ 
ferential inclusions in admit an exhaustive set of zero-separating curves, 
we obtain a proof of the general three-dimensional global attractor conjec¬ 
ture. This was first proved in [Tj, but the construction there is significantly 
more complex. On the other hand, no proof of the general four-dimensional 
global attractor conjecture has been found before, since for that we would 
need an exhaustive set of zero-separating surfaces for three-dimensional toric 
differential inclusions in This is the topic of the next section. 

4.3 Three-dimensional toric differential inclusions 

Consider a toric differential inclusion denoted Tjj in Then, without loss 
of generality we can assume that Ts is generated by a set of planes Tfs such 
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that the set of lines £3 that are orthogonal on the planes in 1 -i^ satishes 
span{C^) = 

We will also assume that for any line I E C^, the lines that can be obtained 
from I by linear transformations given by permutations of the coordinate axes 
are also contained in £3 (there are 3! = 6 such lines). In other words, we 
assume that £3 is symmetric with respect to linear transformations given by 
permutations of the coordinate axes. 

Moreover, we will also assume that the three coordinate axes are con¬ 
tained in £3, and the lines of the form {x = 0,y = z}, {y = 0, z = x}, {z = 
0 ,y = x}, are also contained in £3. 

It is difficult to represent the toric differential inclusion Ts geometrically 
in since we would have to draw many planes and all their intersections 
along lines. To simplify our figures, we will represent Ts geometrically in 
as described below. 

In general, it is easier to describe the associated polar differential inclu¬ 
sion, polar(T 3 ), instead of % itself. Since the connection between % and 
polar(T 3 ) is a simple logarithmic diffeomorphism ([^ we will often describe 
Ts using polar{T 3 ), and will refer to the domain of polar{T 3 ) as logarithmic 
space. 

Since the information about Ts is contained in the set ^{3 of planes that 
bound the regions on which polar{T 3 ) is constant, we describe a way to 
visualize these planes in Note also that only the planes that intersect 
the negative orthant are of interest for the construction of a zero-separating 
surfac^^ Moreover, all these planes contain the origin, so they are in one- 
to-one correspondence with lines in the plane X + Y + Z = —1 that intersect 
the triangle {X + Y + Z = —1} n 

Moreover, due to the permutation symmetries mentioned above, it is 
sufficient if we construct a zero-separating surface only in one of the 6 regions 
given by specifying inequalities between x, y and z, for example in the region 
{x < y < z} n which is shaded in Fig. i K we map this region to 
logarithmic space we obtain the region {X <Y < Z} C which, when 
intersected with the triangle {X + Y + Z = —l}n gives us a smaller 
triangle, similar to the one shaded in Fig. (except that it is in the negative 
orthant). 

Further, we map this smaller triangle (given by {X <Y< Z}r\{X + Y + 

similarly, in the previous section, only the lines that intersected the 3rd quadrant 
were relevant to our construction of a zero-separating curve. 
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Figure 3: We partition the domain into six regions that are mapped into 
each-other by the six symmetries given by permutations of the axes. The 
shaded triangle points out the region IZ^ = {{x,y,z) G | x <y < z}. 

Z = —1} nMl) to the plane {Z = 1} by dividin^^ all the coordinates by Z. 
This transformation can be regarded geometrically as follows: we consider the 
line through a point (X, Y, Z) and the origin; then we map this point to the 
intersection between this line and the plane {Z = 1}. In particular, it follows 
that line segments inside the domain {X <Y < Z}r\{X -\-Y -\-Z = — l}nM^ 
are mapped to line segments inside the plane {Z = 1}; moreover, the image 

we can regard the new coordinates in the plane {Z = 1} are projective coordinates. 
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of this mapping is {X > y > Z} fl {Z = 1} C and can also be written 
as {(X, y, 1) e I X > y > i}. 

Note now that each plane that intersects the region {X < y < Z} of the 
negative orthant can be represented uniquely by a line segment within the 
triangle {X < y < Z} fl {X + y + Z = — l}n (such that the endpoints 
of the line segment are on the edges of the triangle), which in turn can be 
represented uniquely by a line segment or a half-line within the “unbounded 
triangle” {(X, y, |X>y>l}. The endpoints of these segments or 

half-lines are on the edges of the unbounded triangle. 

Moreover, we can map the set {(X,y, 1) G | X > y > 1} to 
by simply projecting on the Xy-plane. Then the toric differential inclusion 
Ts can be represented by a set of line segments and half-lines in the region 
{(X, y) e I X > y > 1}, as shown in Fig. |^a). (We do not need to 
consider all 6 such regions since we have assumed that £3 is symmetric with 
respect to linear transformations given by permutations of the coordinate 
axes.) 

Note that, without loss of generality, we can assume that Ts also contains 
additional horizontal half-lines through all the intersection points of line seg¬ 
ments and half-lines in Fig. |^b), and additional non-horizontal line segments 
connecting these intersection points, such that all the regions bounded by line 
segments are simplicial (i.e., triangular), as shown in Fig. |^c). In particular, 
this means that we now no longer regard Tj as a hyperplane-generated toric 
differential inclusion, but it is now a general toric differential inclusion. 

We consider now the set 1^2 = GM^|X>y>l} and its 

simplicial partition with line segments and half-lines described above and in 
Fig. 1^ Consider also the set TZ^ = {{x,y,z) G \ x < y < z}. We will 
construct a function ip : V 2 ^ TZs whose graph (together with its symmetric 
images with respect to permutations of axes) will form a zero-separating 
surface for Tjj. 

Consider a faithful zero-separating curve (obtained as described in the 
previous section) along the half line {(X, y) G | X = Y > 1} in £> 2 ; 
which is the left-side boundary of £>2 in Fig- Consider the partition of the 
domain 1^2 shown in Fig. This partition was obtained as follows. 

The shaded stripes start on the left side at the points given by the zero- 
separating curve mentioned above, and go approximately horizontally across, 
as shown. We will see later how exactly we choose the upper and lower bound¬ 
aries of the shaded stripe regions. Along each shaded stripe region there are 
vertical line segments in between each pair of consecutive simplicial parti- 
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Figure 4: We partition the domain V 2 into triangles, in several stages: (a) To 
begin with, the domain V 2 is crossed by lines that represent the planes of the 
hyperplane-generated polar{%.)■ Then, in (b) we draw additional horizontal 
lines through all the intersection points of lines in (a). Finally, in (c) we 
draw additional non-horizontal line segments to subdivide each polygon in 
(b) into triangles, without adding additional intersection points. 
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Figure 5: This diagram represents a further partition 0 /P 2 into regions. Our 
construction of the map ip (whose graph is a zero-separating surface for %) 
uses this diagram to proceed from one region to the next, as explained below. 
Regions are denoted Sij, and the curves Ci and Rij represent boundaries 
between regions. Each region Sij contains exactly one point where two or 
more solid lines intersect, denoted Pij. 


tion vertices, and these vertical line segments connect a point on the lower 
boundary of the shaded stripe with another point on the upper boundary 
of the shaded stripe, as shown in Fig. Also, there is exactly one dotted 
line segment within each region of the simplicial partition, and it connects 
an upper boundary point on the lower shaded stripe with a lower boundary 
point on the upper shaded stripe, as shown in Fig. 

We now describe in more detail how this partition is generated, and how 
it allows us to construct the function (p : V 2 ^ TZs = {{x,y,z) G | x < 
y < z} mentioned above. 

Note that the union of curves {(f^, , f) | f > 0} forms a disjoint 

(X,Y)£V2 

partition of the set near the origin, and any zero-separating surface close 
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enough to the origin must contain a point on each one of the curves in the 
union above. 

We will think of each point (X, Y) G T >2 as associated to the curve 
\ t > 0}, and we will dehne the function ip such that ip{X,Y) 
lies on this curve. In other words, consider the function -0 : TZ^ fl —)■ P 2 

such that for each point (x, y, z) G 71^ at distance at most £ to the origin we 
have that 'ip{x,y,z) = {X,Y) G "^2 such that there exists some t > 0 with 
{x,y,z) = (Note that actually t = z.) Then, we will dehne the 

function ip such that 'ip o p> ■. T >2 ^ T >2 is the identity function. 

This means that if we specify a point (xq, 1 / 0 , 2 ^ 0 ) ^ which is on the 
graph of the function (p, then we have also implicitly specihed exactly for 
which point (Xq, To) G ”^2 we have </5(Xo, Xq) = (a^o, UoiZo). So, to dehne the 
function ip, it is enough to specify the surface in TZ^ which is the graph 
of (p. 

As we mentioned above, we start by dehning the graph of (p along the 
left-side boundary of 'D 2 , i.e., along the half-line {(X, X) G | X = X > 1}, 
and we do so by using a faithful two-dimensional zero-separating surface (i.e., 
curve) along this half-line, constructed as described in the previous section. 
This gives us a left-side boundary for the surface S^p (if we think of the surface 

as lying in Fig. where the plane {x = y} is to the left of the region 7 ^ 3 ). 

In particular, we choose the flat regions of the left-side boundary of 
to be short enough, so the shaded regions are thin enough to allow the ver¬ 
tical lines shown in each shaded region to cross that shaded region without 
intersecting other lines (we may have to move the construction closer to the 
origin to make this possible )p^ 

We now explain how we dehne (p on the various patches of 1)2 shown in 
Fig. and, equivalently, we construct patches of the surface in TZ^, such 
that when we map each such patch to T >2 via "0, we obtain one of the patches 
of 7)2 shown in Fig. 

The left-bottom patch of S^p is hat, given by a plane whose normal vector 
is (Pii, 1) = (1,1,1), where the point Pn is the corner of the domain V 2 - 
Note that this hat patch 5ii is compatible with the zero-separating curve on 
the left-side boundary of S^. 

^®We plan to construct a zero-separating surface based on the “blueprint” shown in 
Fig.0 The blueprint comes with a specified width for the shaded horizontal stripes, to 
accommodate other features, as we will see later. We don’t need to change this blueprint 
to accommodate a larger value of <5; instead, we can compensate for the larger 5 just by 
moving the whole construction of the zero-separating surface closer to the origin. 
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The boundary of this patch is made up of line segments, and we describe 
them as follows. Note how in Fig. the boundary of the image of this 
patch in V 2 consists of three curves, denoted Ci, (72,6*3. Each one of these 
curves intersects exactly one solid line segment, and recall that the solid 
line segments in Fig. [^represent the separating surface^^oi T 3 . To simply 
the terminology we will say that to each black line segment in Fig. there 
corresponds an attracting directiotf^ of Ts. 



Figure 6: Zoom in on the left-side corner of the diagram in Fig. 

Then, we choose Ci such that the boundary ip{Ci) of the flat patch 
5ii is a line segment in the attracting direction of the black line segment 
that Cl intersects, and we choose C 2 and similarly. Since each one of 
ip{Ci),ip{C 2 ),'p{Cs) is orthogonal to a plane that contains the line corre¬ 
sponding to the point (Pn,!), it follows that they do indeed belong to a 
common plane, which is orthogonal to the vector (Pn, 1). 

Note that, since (^((Fs) is a line segment in the attracting direction of 
a horizontal segment in Fig. it follows that the a:-coordinates of points 
along (^((Fa) must be constant; therefore the X-coordinates of points along 
(F3 must also be constant, because (F3 = fj{ip{C^)). Therefore, (F3 is flat and 
orthogonal on the horizontal black line segment it intersects. 

^®the separating surfaces of Ta are exponential images of the hyperplanes that generate 
polar{%). 

As we discussed before, the attracting directions are easier to describe for polar{Tf)', 
but recall that if we take into account the correspondence between their separating sur¬ 
faces, then Ts and polarifC) have the same attracting directions. 
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On the other hand, ^{ 02 ) is a line segment in the attracting direction of 
a non-horizontal segment in Fig. so the x-coordinates along 93 ( 02 ) are not 
constant. This means that O 2 is not perfectly horizontal, although it may 
appear so in Fig. and Fig. But, when we map 93 ( 02 ) through ip to obtain 
O 2 , the image must be a curve with almost horizontal tangent line. We can 
see this by checking that, if we £x some values X > Y > 1 and 0 < f -C 1 
and we calculate the Jacobian of ip at the point {t^,t^,t), we obtain the 
matrix 


log^ 

t 


0 -X \ 

0 0-’^ -Y ) • 


( 10 ) 


Note now that, given the assumptions above we have t^~^ 3> X and t^~^ 3> 
y, and also t^~^ t^~'^. This explains why O 2 is approximately horizontal 

everywhere. The case of Oi is similar, except that Oi reaches the line {X = 
y} in V 2 , so Cl is approximately horizontal far enough from that line, while 
close enough to the line {X = y} the curve Ci becomes approximately 
orthogonal to it. 

Let us look now at the second shaded patch 5 i 2 on the bottom shaded 
stripe, such that ip{Si 2 ) is bounded by the curves € 3,04 and C 5 in Fig. 
This will also be a flat patch, given by a plane whose normal vector is (P 12 ,1), 
where P 12 is the intersection point of solid lines within the patch 1 S 12 . 

Note hrst that the patch 5 i 2 is compatible with the patch iSn we con¬ 
structed before, because they meet along the line segment 93 ( 173 ), and, since 
93 ( 173 ) is orthogonal to the plane in determined by the horizontal black 
line it intersects, it follows that it is also orthogonal to any line within this 
plane, and in particular it is orthogonal to the vectors (Pn, 1 ) and (P 12 , !)■ 
Therefore the patch 5 i 2 is compatible with the patch ^n, i.e., they intersect 
along the line segment 93 ( 173 ) whose direction is indeed the intersection of 
the two planes that give and 5 i 2 . 

The construction of C 4 is similar to the construction of C 2 , and the con¬ 
struction of C 5 is similar to the construction of C 3 . 

After this, the construction of the other patches Sij in the on the bottom 
shaded stripe of Fig. [^proceeds similarly. We will now discuss the construc¬ 
tion of the patches along the bottom white stripe of Fig. Let us start by 
looking at the patch denoted 521 . 

As before, we start with a zero-separating curve along the left-side bound¬ 
ary of this patch. Starting from each point of this curve, we draw a line 
segment in the attracting direction of the left-side boundary line, and con- 
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tinue it to the right until it reaches the surface represented by the dotted line 
segment R21 in Fig. We obtain a curve along the surface represented by 
and then, starting from each point of this curve, we draw a line segment 
in the attracting direction of the next black boundary line, and continue it to 
the right until it reaches the surface represented by the dotted line segment 
R22 in Fig. This will give us a curve along the surface represented by i?22, 
and we continue the construction in the same way until we hll in this bottom 
white stripe of Fig. [^completely. 

Next, we construct the zero-separating surface ^(T>2) along the second 
shaded stripe of Fig. [^ then the second white stripe, and so on. Along each 
stripe the construction starts on the left-hand side of the domain V2 and, 
one by one, hlls in each region in the partition shown in Fig. 

We need now to check that ^{^2) (together with its 6 symmetric images) 
is indeed a zero-separating surface for Ta. For this it is enough to see that, 
within each region in the partition shown in Fig. [^ the surface ‘^{^2) is a 
support surface of the cone of Ts in that region. 

First we check that ^{^2) satishes the conditions necessary for being 
a zero-separating surface for Ts along the shaded region patches. Indeed, 
all the surface patches of (p(T>2) in the shaded regions of V2 are flat plane 
patches, and for each shaded patch that plane is exactly the cone of Ts at 
the black point that lies at the center of that patch, so it follows that ^{^2) 
is a support surface of the cone of % in all shaded regions in Fig. 

Consider now the hrst white patch we constructed, denoted 52i. We will 
argue that the tangent plane at any point of this patch is a support plane 
of the cone of Ts in that leftmost lower triangle of the partition in Fig. [^ 
Recall that the left half of ^21 (up to the surface represented by R21) was 
constructed as a ruled surface that started from the left-side boundary curve 
Co, which was part of a faithful zero-separating curve along the line {X = R}. 
The direction of the straight line segments along this ruled surface was the 
attracting direction of the surface represented by T?2i- Then the tangent 
plane to this ruled surface must contain this constant attracting direction. 

If we imagine the construction above for the case where the left-side 
boundary curve Cq is replaced by a short line segment /q within the plane 
{x = y}, then the ruled surface is a plane patch tto. Then note that the 
normal to tto is the same as the normal to lo within the plane {x = y}, 
because tto is orthogonal to {x = ?/}|^ Therefore, by approximating Cq with 

general, if we have a codimension-two linear subspace I contained in a hyperplane 
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polygonal lines and then taking the limit, we obtain that the normal to 52 i 
at a point (p{P) is approximately the same as the normal to Cq within the 
plane {x = y}, calculated at a point Pq which is obtained by drawing a 
horizontal line from P and intersecting it with the line {X = F} in Fig. 
(here we are using the fact that the image via 'ip oi a line segment that passes 
through (p{P) and is orthogonal on {x = y} is approximately the horizontal 
line from P to the line {X = y}). Then we obtain that the surface patch 52 i 
does satisfy the condition that its outer normal vector belongs to the cone 
corresponding to the region V’(‘ 52 i) in T>2, so 1S21 is a zero-separating surface 
patcl^ for Ts- 

Consider now the second white patch we constructed, denoted ^22 in 
Fig. [ 5 j The surface represented by R21 is asymptotically close to the plane 
{x = 0 } in (i.e., if the neighborhood of the origin in Vq C is small 
enough then this surface is arbitrarily close even in relative terms) 

Then, consider a polygonal line approximation of the right-side boundary 
curve of ^21, such that each line segment in this polygonal line is approxi¬ 
mately parallel to the plane {x = 0 }. Denote one such polygonal line segment 
by /i. If we construct a plane starting from li and going in the attracting 
direction of the surface represented by the solid line in region S22 in Fig. 
then the normal to this plane will have to lie on the plane represented by 
the solid line and also will have to be contained in the dihedral angle rep¬ 
resented by the white stripe (because it has to be orthogonal both to the 

%, and we build a hyperplane tt by extending from I in a direction orthogonal to R, then 
the normal line to I within R is the same as the normal line of the hyperplane tt. (This 
will be useful for higher dimensional constructions similar to the one above, with I inside 
the hyperplane R = {x = y} in M".) 

i.e., it is a surface patch that has the local properties of a zero-separating surface for 

%. 

Consider a surface z = x‘^y^ that contains the curve ,t)\t > 0 } for 

some 7i > 72 > 1 - The normal to the surface is given by grad{x°‘y^ — z) = 
{ax°‘~^y^, Px°‘y^~^ ,—l), so the normal at a point Rf) is given by the vector 

,-l) = - 1 )? because 071-1-/372 = 1 . If 

<0 is small enough then 1 - Therefore, if o ^ 0 , the normal vector 

,— 1 ) is approximately along the x-axis. 

This calculation does not include z-independent surfaces of this kind, i.e., surfaces of the 
form 1 = x°‘y^ such that 071 -|- /372 = 0 for some 71 > 72 > 1 - Note that the normal vector 
(atQ~"’'bsimply replaced by ,Q) in this case. Alternatively, 

we can analyze all cases at once by looking at the surface 1 = x°‘y^z'^ instead of the surface 
z = then the normal at the point Rq) becomes and we 

still have 
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attracting direction of the surface represented by the line Cg, and to the line 
segment by li we started with, and because the original zero-separating curve 
that gave us the left-side initial condition was faithful). Therefore the second 
white patch we construct (denoted by S22 in Fig. is also a zero-separating 
surface patch for T2. 

We can now continue towards the right in the lowest white stripe in Fig. 
and conclude that we have constructed a zero-separating surface patch for 
T2 along this whole white stripe in V2. Then we do the same for the next 
higher shaded stripe, and then for the next higher white stripe, and so on. 
In conclusion, it follows that (p(T>2) is a zero-separating surface for 72 . 


Note that the zero-separating surface (p(T>2) we have constructed above is 
not faithful. We can obtain a faithful zero-separating surface f)(T>2) by using 
a similar construction, but based on a partition of V2 like the one shown in 


Fig. 10 



Figure 7 : A different partition of V2, that can be used to build a faithful 
zero-separating surface for Tg. 


Note that each triangular patch in Fig. 10 is now divided into three sub¬ 
patches (and not just two as before). Of these three sub-patches, we dehne 
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(p along the left-side and right-side patches as before, but for the middle 
sub-patch we proceed differently, to ensure that we obtain a faithful zero- 
separating surface. 

For each middle sub-patch within a shaded stripe we dehne (p to be not linear 
along the horizontal direction (as before) but to vary along this direction in 
such a way that the normal to the graph of at a point (p{P) is very close 
to the projection of P on the black horizontal line along the patch. (This 
property can be ensured by first dehning (p along the black line, and then 
extending the graph of (p to the middle sub-patch via line segments that are 
orthogonal to the black line.) 

Also, for each middle sub-patch within a white stripe we define (p such that 
the normal to the graph of at a point (p{P) is very close to P itself. (One 
way to accomplish this is formulate this condition as a hrst order PDF with 
left-side boundary condition given by the already constructed graph of (p on 
the left side of the middle sub-patch.) 


4.4 Four-dimensional toric differential inclusions 

Consider a toric differential inclusion denoted 71 in Without loss of 
generality (analogous with the previous section) we can assume that 71 is 
generated by a set of hyperplanes 7^4 with a set of orthogonal lines £4 such 
that span{Ci) = 

Also analogous with the previous section, we will assume that for any line 
I & C4, the lines that can be obtained from / by linear transformations given 
by permutations of the coordinate axes are also contained in £4 (there are 
4 ! = 24 such lines). In other words, we assume that £4 is symmetric with 
respect to linear transformations given by permutations of the coordinate 
axes. 

Moreover, we will also assume that the three coordinate axes are con¬ 
tained in £4, and the lines of the form {x = y, z = 0,w = 0}, {x = z,y = 
0 ,w = 0}, {y = z,x = 0,w = 0}, {x = w,y = 0, z = 0}, {y = w,x = 0, z = 
0 }, {z = w,x = Q,y = Q} are also contained in £4. 

It is difficult to represent the toric differential inclusion 71 geometrically 
in in a useful way. To be able to draw some useful hgures, we will represent 
key parts of 71 geometrically in as described below. 

All the information about 71 is contained in the set of hyperplanes that 
bound the cones of its associated polar differential inclusion, po/ar( 71 ). Only 



the hyperplanes that intersect the negative orthant are of interest for the con¬ 
struction of a zero-separating surface. Moreover, all these hyperplanes con¬ 
tain the origin, so they are in one-to-one correspondence with two-dimensional 
planes in the hyperplane X + Y + Z + W = —1 that intersect the tetrahedron 

{x + r + z + hh = -i}nMi. 



Figure 8: Like in the previous section, we use the symmetries ofT^ to divide 
the problem into 41=24 smaller problems. The domain is obtained by 
mapping into one of the 24 tetrahedra, such as the one with thick edges 
shown above. 

Moreover, due to the permutation symmetries mentioned above, it is 
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Figure 9: The construction of the function ip in the previous section has 
proceeded one step at a time, first the bottom shaded layer of V 2 , then the 
bottom white layer, and so on. The analogous construction in this section 
features 3d “tunnels” instead of 2d layers, as shown above. In the figure above 
we should think of the configuration of solid lines as being all in the same 
plane, on the left-side facet ofV^, and representing the domain of a left-side 
boundary condition for ip. Then, each solid line segment of the left-side facet 
ofV^ determines a (flat) surface patch inside V^, obtained by drawing lines 
in the X-direction that start on that line segment. The dotted lines shown 
above are lines of this kind that happen to pass through intersection points 
of two or more segments on the leftside facet ofV^. These dotted lines and 
the flat surface patches that connect them determine the boundaries of the 
3d tunnels we mentioned above. There is additional subdivision within these 


patches, as explained in Fig. 10 
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sufficient if we construct a zero-separating surface only in one of the 24 regions 
given by specifying inequalities between x, y, z and w, for example in the 
region {x < y < z < tc} If we map this region to logarithmic space we 

obtain the region {X < Y < Z < W} <Z which, when intersected with 
the tetrahedron {X-|-y-|-Z-|-hh = gives us a smaller tetrahedron, 

similar to the situation shown in Fig. 

Further, we map this smaller tetrahedron (given by {X < Y < Z < 
IF} n {X + Y + Z + W = —1} n Mi) to the plane {IF = 1} by dividing all 
the coordinates by IF. This transformation can be regarded geometrically 
as follows: we consider the line through a point (X, Y, Z,W) and the origin; 
then we map this point to intersection between this line and the hyperplane 
(IF = 1}. In particular, it follows that line segments inside the domain 
(X < Y <Z<IF}n (X -I- F -I- Z -I- IF = —1} n Mi are mapped to line 
segments inside the plane (IF = 1}; moreover, the image of this mapping is 
(X >Y>Z> IF}n{IF = 1} C M{_, and can also be written {(X, F, F, 1) G 
M^ I X > F > IF > 1}. 

Therefore each hyperplane that intersects the region (X < Y < Z < 
IF} of the negative orthant can be represented uniquely by the intersection 
between a two-dimensional plane and the tetrahedron (X < Y < Z < 
IF}n{X-|-F-l-F-l-IF = — l}nMi, which in turn can be represented uniquely 
by the intersection between a two-dimensional plane and the “unbounded 
tetrahedron” {(X, F, F, 1) G M"^ | X > F > F > 1}. Moreover, we can map 
the set {(X, F, F, 1) G M^ | X > F > F > 1} to M^ by simply projecting on 
the XFF-space. The result is shown in Fig. [^and Fig. 

Then the toric differential inclusion 74 can be represented by a such a 
set of planes in the region {(X, F, F) G M^ | X > F > F > 1}. (We 
don’t need to consider all 24 such regions since we have assumed that £4 is 
symmetric with respect to linear transformations given by permutations of 
the coordinate axes.) 

Note that, without loss of generality, we can assume that 71 also contains 
additional horizontal “unbounded triangles” {(X, F, F) G M^ | X > F > F > 
1, F = constant} through all the intersection points of planes in Fig. and 
additional non-horizontal triangles connecting these intersection points, such 
that all the regions bounded by line segments are simplicial (i.e., tetrahedral), 
as illustrated in Fig. 10 (Like in the previous section, this means that 71 


is no longer a hyperplane-generated toric differential inclusion, but is now a 
general toric differential inclusion.) 


31 




Figure 10: A partition of V^, analogous to the partition of V 2 from in the 
previous section, uses tetrahedra instead of triangles. Here we show only 
two representative “tunnels” in this partition, and the six different types of 
tetrahedra that may occur, depending on the positions of the vertices of the 
front face and the back face. Note that in Figl^c) we also had two types 
of triangles (with the horizontal edge on the bottom or on the top). This 
difference is reflected here in the difference between the top three and the 
bottom three tetrahedra. 

In order to describe a method that allows us to construct a zero-separating 
surface for 71, let us revisit the construction of a zero-separating surface for 
71 in the previous section, to emphasize some key ideas. 

If we look again at Fig. we can interpret it as follows. On the left¬ 
side boundary of the domain V 2 we had two types of curve patches, i.e., the 
Id-flat patche^^ and the Od-flat patched 

Then, if the Id-flat patches are short enough, and if the whole left-side 

By Id-flat patches we mean the parts of the left-side boundary curve near some 
vertices, which give rise to the shaded layers in the partition of T>2 shown in Fig. The 
proximity of the vertices impose that these curve patches are straight line segments. 

By Od-flat patches we mean the parts of the left-side boundary curve that are not close 
to vertices, and give rise to the white layers in the partition of 'D2 shown in Fig. Some 
of these patches may actually be flat (i.e., straight line segments), but there is nothing in 
our construction that requires them to be flat. 
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boundary curve is faithful, then the construction described in the previous 
section is possible. 

We can think of this construction systematically as follows. For each Id- 
flat patch we have constructed a polygonal line of Id-flat patche^^ which is 
the image through the function ip of the corresponding shaded layer in the 
partition of V 2 shown in Fig. Similarly, for each Od-flat patch we have 
constructed a polygonal line of Od-flat patches. For both types of polygo¬ 
nal lines of patches the construction is similar: the patches are continued 
from left to right in the attracting direction of a surface represented by a 
solid line segment, and they “turn” in a different direction along the surface 
represented by a dotted line segment. 

A key property that makes this construction work is the fact that the 
two different types of polygonal lines of patches match perfectly along their 
common boundaries. On the other hand, we can see in Fig. that the 
polygonal line of Od-flat patches “turns” (i.e., encounters dotted lines) more 
times than the polygonal line of Id-flat patches. The reason their boundaries 
still match is the following: even though the polygonal line of Id-flat patches 
does not turn as many times (e.g., does not turn in the neighborhoods of 
points where solid lines intersect), its boundary curves do turri^^ 

Note that all we needed in the previous section in order to construct the 
polygonal lines of patches described above were the two different kinds of 
dotted lines. One way to think about our choice of a conhguration of dotted 
lines as shown in Fig. is the following. We can hrst choose “midpoints” for 
all the horizontal solid line segments, and decide that the thick dotted lines 
will pass through these midpoints. Once we have the thick dotted lines, we 
can also construct the intersections between them and the thin horizontal 
dotted lines. These intersection points are important, since the only other 
remaining dotted lines (the non-horizontal ones) pass through these points. 
We can construct these lines for example by connecting the intersection points 
with the triangle vertex that lies on the opposite side. 


We call it a polygonal line of Id-flat patches because we construct it in a way that 
is similar to a polygonal line, except that a polygonal line starts with a point, and then 
advances along some straight lines, while a “polygonal line of Id-flat patches” starts with 
a Id-flat curve Ci, then builds a surface patch by advancing along straight parallel lines 
until we reach another curve C2 at the intersection with some surface, then builds another 
surface patch by advancing along another set of straight parallel lines, and so on. 

and even though the boundary curves turn, they still remain part of the same Id-flat 
patch. 
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We now construct a zero-separating surface for 71 based on the idea of 
polygonal lines of patches. We need to describe “dotted surfaces” (i.e., sur¬ 
faces that sub-partition D 3 , and play the role that the dotted lines have 
played for sub-partitioning P 2 in Fig- for the tetrahedral partition of the 
three-dimensional domain P 3 illustrated in Fig. |^and Fig. 

We construct the dotted surfaces as follows. First, like in the 3d case, 
we choose midpoint^^ for all the horizontal line segments in the tetrahedral 
subdivision of For example, in Fig. 10 (where we show just two repre¬ 
sentative tunnels within V 3 ) there will be one such midpoint on each line 
segment along the six horizontal half-lines shown there. Note that, in order 
to make this hgure more clear, we have omitted many other tetrahedra that 
are also inside these same two tunnels (since, according to our construction, 
7^3 is partitioned by tunnels, and each tunnel is partitioned by tetrahedra). 

Then, once we have the midpoints, we also construct mid-curvei^ along 
the X-direction tetrahedral face^^ Each mid-curve connects a point on the 
boundary of its adjacent midpoint turning pafc/Q with the opposite vertex 
of the X-direction tetrahedral face that contains the mid-curve. 

Once we have specified the midpoints and the mid-curves as described 
above, there is a unique way to build the polygonal lines of 2 d-flat patches 
and the polygonal lines of Id-flat patches. 

Finally, we construct a mid-surf ac^^wiihin each tetrahedron, and it gives 
us a unique way to build the polygonal lines of Od-flat patches, and therefore 
hnish the construction of a zero-separating surface for the toric differential 
inclusion 71. 


^®we need these midpoints to be able to specify the turning points of the polygonal lines 
of 2 d-flat patches that we will construct later. 

the mid-curves will specify the turning points of the polygonal lines of 1 d-flat patches 
The X-direction tetrahedral faces are faces that contain the direction of the X-axis. 
Each tetrahedron contains exactly two such faces, and they are adjacent to a common 
horizontal edge, so they both have an adjacent midpoint that we constructed previously. 
The other two faces are called back face and front face, as seen in the direction of the 
construction, i.e, in the X-direction. The back face and the front face intersect in an edge 
which is called the hinge edge of the tetrahedron. 

a midpoint turning patch is the vertical patch that passes through a midpoint and 
along which the polygonal line of 2 d-flat patches turns. Each X-direction tetrahedral face 
intersects exactly one midpoint turning patch. 

^^the mid-surface of a tetrahedron is a surface patch that lies approximately half-way 
between the back face and the front face of the tetrahedron, and it connects the hinge 
edge with the boundaries of the polygonal lines of Id-flat patches and 2 d-flat patches that 
intersect that tetrahedron. 
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As a final remark, let us also notice that there is another way to describe 
the above construction, as follows. For each tetrahedron in the subdivision 
of Vs, we choose a midpoint of its advance and then construct two 

“dotted lines” that start at this midpoint and go into the two faces of the 
tetrahedron that contain the advance edge Then the “dotted surface” 
of the tetrahedron is a triangular surface ^ that partitions the tetrahedron 
into two parts^ Two edges of the dotted surface are the dotted lines we 


mentioned above, and the third edge is the hinge edge of the tetrahedron. The 
dotted surface hlls in the space in between its three edges according to the 
following rules: {i) the midpoint must be contained in the midpoint turning 
patch of the polygonal lines of 2d-flat patches, and the part of the midpoint 
turning patch that is contained in this tetrahedron must be contained in the 
dotted surface, {ii) the mid-curve must be contained in the mid-curve turning 
patch of the polygonal lines of Id-flat patches, and the part of the mid-curve 
turning patch that is contained in this tetrahedron must be contained in the 
dotted surface, and (Hi) the dotted surface is “vertical enough”, in the sense 
that approximately horizontal line segments intersect it in at most one point. 

Condition (i) is trivial to enforce. To enforce condition (ii) we have to 
make sure that we do not allow the two mid-curve turning patches to intersect 
near the hinge edge. Finally, notice that (by moving the construction of the 
zero-separating surface close enough to the origin) we can assume that the 
turning patches as as thin as we wish, so the fitted surface is as close to the 
half-distance triangle between the back face and the front face. Therefore 
condition (Hi) can also be satished. 


Our next focus will be to construct a faithful zero-separating surfac^^for 
the 4d toric differential inclusion T 4 . 

^^the advance edge of a tetrahedron is its only edge in the horizontal direction (i.e., in 
the direction of the X-axis). 

^"‘these two faces are the horizontal faces of the tetrahedron. 

^^it is triangle-shaped, but it is not flat. 

Recall that the dotted surface has to serve as “turning surface” for the three kinds of 
polygonal lines of patches. The dotted surface has to partition the tetrahedron into two 
parts, one containing the back face, and one containing the front face. Then we have a 
“back half’ and a “front half’ of the tetrahedron, and in the back half the polygonal lines 
of patches advance in the attracting direction of the back face, and in the front half the 
polygonal lines of patches advance in the attracting direction of the front face. 

because it will be the starting point for a construction of zero-separating surface for 
a 5 d toric differential inclusion T5, in the next section. 
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Recall that we have constructed a faithful zero-separating surface for the 
3d toric differential inclusion Ta in the previous section, essentially by build¬ 
ing a slightly smaller triangle inside each triangle of the subdivision of 1 ^ 2 , 
in a special way that makes all of the vertices of these smaller triangles align 
with each-other along turning patches and along some attracting directions 
(see Fig. 10). 

Surprisingly, although this simple approach worked well in it fails to 
worl^ in Therefore, we need a new way of constructing a faithful zero- 
separating surface for 71, which will also generalize to higher dimensions. 

A method that accomplishes this is the following. After we have subdi¬ 
vided the domain T>^ of 71 into tetrahedra, we further subdivide each tetrahe¬ 
dron by choosing several points along its advance edge and connecting them 
with its hinge edge, and making sure that the hrst such point is very close 
to the back face, the last such point is very close to the front face, and there 
is at least one more point between them. 

Denote this new polygonal fan by 7^1, and its associated toric differential 
inclusion by Tf, and consider a zero-separating surface for Tf constructed 
close enough to the origin. Then, for each tetrahedron of 71, the normals 
to this zero-separating surface for 72 at points that are not in any uncer¬ 
tainty region of 71 belongs to the interior of the attracting cone of that 
tetrahedron of 71. Therefore, this zero-separating surface for Tf is a faithful 
zero-separating surface for 71. 


4.5 General (n-dimensional) toric differential inclusions 

The constructions we have described in the previous sections can be carried 
on analogously in higher dimensions. To build a zero-separating surface for an 
n-dimensional toric differential inclusion 71, we first subdivide the simplicial 
domain 

{Xi + A2 + ...+X„ = -l}nM2 

into n! subdomains given by specific orderings of the logarithmic coordinates 
Xi, ...,Xn- Then we take the representative subdomain 

{Xi < X 2 < ... < xj n {Xi + X 2 +... + X, = -1} n Ml 

and map it to the unbounded simplicial domain C M^”^. 

It fails because the system of restrictions given by these alignments is overdetermined 
in 
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Then we subdivide T’n-i into (n —l)-dimensional simplices (similar to our 
decomposition of P 2 into triangles and of V 3 into tetrahedra) and assume 
that we have already constructed a “left-side boundary condition” for a zero 
separating surface on 

The simplices of can be organized along “tunnels” that start at 

simplices on the left-side boundary of and go in the Xi-direction. This 
provides a lexicographic-type ordering of all the simplices of 'Dn-i, based on 
the ordering of simplices of T>n -2 and the order of increasing Xi values in 
simplices within the same tunnel in T>n-i- 

Then, we build a zero-separating surface for Tn by constructing one piece 
of the surface for each simplex of in the order described above. The 

construction is based on polygonal lines of patches, like in the previous sec¬ 
tion. The turning sets for the polygonal lines of patches are given by a 
“dotted surface” constructed like in the previous section: for each simplex 
of Vn-i, we hrst choose a midpoint of its advance edge, then we construct 
boundary patches of the dotted surface on all the facet^of the simplex that 
contain the advance edge, and then connect these boundary patches with the 
codimension-two face that is opposite the advance edge. 

Finally, we can also obtain faithful zero-separating surfaces for Tn, by 
subdividing each simplex of T>n-i into several simplices, given by a choice 
of at least three points along the advance edge of the simplex, as explained 
in the previous section. Faithful zero-separating surfaces for n-dimensional 
toric differential inclusions serve as left-side boundary conditions for zero- 
separating surfaces for [n + l)-dimensional toric differential inclusions, and 
so on. 


5 Proof of the global attractor conjecture 

Consider a toric dynamical systemp] in M”, and £x some 0 < £ << 1. 
From the work of Horn and Jackson [Hj we know that has a globally 

If we partition into half-lines in the direction of the Xi-axis, then the “left-side 

boundary” of is the set of boundary points of these half-lines. The left-side boundary 

of Vn-i is essentially a copy of 'Dn-2, and the zero-separating surface constructed along it 
is based on previously constructed faithful zero-separating surface for a (n — 1)-dimensional 
toric differential inclusion. 

this relies on the construction of the dotted surface in the (n — l)-dimensional case 
®^i.e., a dynamical system of the form Q that has a positive equilibrium that satisfies 
the vertex-balance identity (p|. 
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defined strict Lyapunov function within any linear invariant subspace. Also, 
we know [ 8 ] that there exist level sets of this Lyapunov function that allow 
us to build an invariant region TZ^ for T„ in such that 

{x G 7?.° 0 ) < £o} = 0 and {x G 7?.° 0 ) > l/^o} = 0 , 


for some £9 > 0, and such that the hypercube "Hg = [e, -]"■ is contained inside 

^n' 

Then we proceed as in the Section 7 of jTj, and build an invariant region 
(Z 77[^ that also contains 77^ and does not have any points at distance 
less than £1 > 0 from the coordinate axes of M”, by using the fact that, in a 
neighborhood of each coordinate axis of , the n-dimensional toric dynam¬ 
ical system restricted to 77° can be regarded as an {n — l)-dimensional 
fc-variable toric dynamical system. This allows us construct 77^ by excluding 
n cylindrical neighborhoods of the coordinate axes from 77°. See also Fig. 7.1 
in d]. 

In the next step we build an invariant region 77^ C 77,'^ that contains "Hg 
and does not have any points at distance less than £2 > 0 from the coordinate 
planes of M", and then we build an invariant region 77° C 77^ that contains 
'He and does not have any points at distance less than £3 > 0 from the 
coordinate 3-spaces of M":, and so on. 


In the enc°^ 
contains 77 


we obtain a bounded invariant region 77” ^ for T„ that 


and does not contain any points at distance less than Sn-i > 0 
from the boundary of M”. Therefore, for each initial condition xq G M” we 
can construct an invariant region 77 ”“^ above such that xq G 77 ”“°. Then 
we can use LaSalle invariance [ElEo] for the global strict Lyapunov function 
mentioned above to conclude that, within any linear invariant subspace, any 
trajectory of T„ converges to its vertex-balanced positive equilibrium point. 
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